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Abstract 

We construct a transformation that transforms perturbative states into states 
that implement Gauss's law for 'pure gluonic' Yang-Mihs theory and QCD. 
The fact that this transformation is not and cannot be unitary has special 
significance. Previous -work has shown that only states that are unitarily 
equivalent to perturbative states necessarily give the same S-matrix elements 
as are obtained with Feynman rules. 
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In earlier work, one of us (KH) quantized Yang-Mills theory in the temporal {Aq = 0) 
gauge and formulated the constraint that implements Gauss's law by selecting an appropriate 
subspace for the dynamical time evolution of state vectors One objective of Ref. was 
to compare the implications of Gauss's law when it is imposed on the non-Abelian Yang- 
Mills (YM) theory and Quantum Chromodynamics (QCD) with its role in QED. Despite 
great similarities between the Abelian and the non-Abelian theories, the inclusion of gauge 
fields in the non-Abelian charge density is responsible for important differences between 
QED and QCD [0,0. Not only is it far more difficult to construct states that implement 
Gauss's law in non-Abelian gauge theories than in QED; it is also much more important to 
use states that implement Gauss's law in evaluating S-matrix elements in YM theory and 
QCD than it is in QED. 

In gauge theories, it is standard practice to use Feynman rules in perturbative calcula- 
tions; these rules implicitly use charged particle states that do not obey Gauss's law. In the 
evaluation of S-matrix elements in QED, perturbative states that do not implement Gauss's 
law may be safely substituted for states that do implement it. This has been shown to be 
due to the fact that unitary transformations suffice to construct the latter states from the 
former But this unitary equivalence does not extend to the non-Abelian gauge theo- 

ries. The validity of perturbative calculations based on Dirac spinor quarks and free gluons 
may therefore require qualifications in YM theory and QCD that are not needed in QED. 
Furthermore, the question has been raised whether the proper implementation of Gauss's 
law in non-Abelian gauge theories might have significant implications for the confinement 
of colored states and the conjectured requirement that only color singlets can be asymptotic 
scattering states [|l|]. 

In this paper, we construct states that obey the non-Abelian Gauss's law in 'pure gluonic' 
YM theory and QCD. Our program is based on the construction of a transformation T 
(which must be non-unitary) that transforms perturbative states \a) — in the first instance 
the perturbative vacuum state |0) — into states that satisfy Gauss's law, and that continue 
to satisfy Gauss's law even after dynamical time evolution. Unlike James and Landshoff, who 
had to require matrix elements of a non-terminating progression of powers of the "Gauss's 
law operator", Jq —diE"-, to vanish in order to obtain states that implement Gauss's law [Q, 
we find that we do not need to apply progressively escalating powers of projection operators 
to achieve our objective. 

As in Refs. [^0], we represent the transverse gauge fields and their adjoint momenta as 
^T.(r)= E ^K(k)e* + af(k)e-^'^-], (1) 

k;s=l,2 V^fc 

and 

^Ur)= E -^em^[a:{k)e-^'^--af{k)e-^^-], (2) 

k;s=l,2 

where a^^h) and a"(k) represent 'standard' creation and annihilation operators, respec- 
tively, for gluons (or — with the Lie group indices removed — photons) of helicity s. The 
longitudinal fields are represented in terms of ghost excitation operators, in the form 
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k 2A;2 

and 

^Ur) = E ^ K(k)e+^''- - a-(k)e-^'^-] . (4) 

The ghost excitation operators obey the commutation rules [aQ(k), a5^(k')] = 5afe5k,k') and 
[a^(k), ag (k')] = 5afe^k,k'; with all other commutators vanishing. We define the "Gauss's law 
operator" ^"(r) 

^"^(r) = 9.n^(r) + Jo"(r) = -^.^^(r) + JS{r) ; (5) 
where Jo{r) = g f°'^'^A'^{r) (r). ^"(r) can also conveniently be represented in the form 

^'^(r) = W [^"(k) + ^^'^^(-k)] 6+^'^- , (6) 
^ k 

whereQ 

n'^{k) = 2kla^Q{k) + JS{k). (7) 

In earlier work [|l],H, it was demonstrated that Gauss's law and the gauge choice, Aq = 0, 
could be implemented by imposing 

Q%k)\u) = (8) 

on a set of states, r2''(k) — and its adjoint r2"*(k) — commute with the Hamiltonians 

for YM theory and QCD, as well as (with the Lie group index a removed) for QED. Gauss's 
law, once imposed by this method, therefore is unaffected by time evolution and remains 
permanently intact. 

In QED, the Noether current commutes with the gauge field, because the photons couple 
to, but do not carry the electric charge. fi(k) therefore commutes with f2*(k), mirroring 
the commutation rule between aQ(k) and aqik). That makes it possible to establish a 
unitary equivalence between the set of states {|z^)} — the solutions of Eq. (|) — and the 
set of states that solve aQ(k)\n) = 0. We are able to exploit this unitary equivalence to 
explicitly construct the states in and to reformulate QED as a theory of charged 

particles that obey Gauss's law — that therefore carry their Coulomb field with them — 
and that interact with each other and with transversely polarized propagating photons . 
In YM theory and QCD, however, the commutation rules among the components of ^"(k) 
are quite different from the commutation rules among the corresponding ^^(k). The latter 
'ghost' annihilation operators commute with each other and with all their conjugate aQ*(k). 
However, the commutation rules for fi"(k) and f2'^*(k) follow the SU{N) Lie algebra [|I|. 



J7"(k) — and rj'**(k) — in this work and in Ref. [||] differ by a normalization factor of 2k^/^. 
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It is therefore impossible to construct a unitary transformation that transforms ^'^(k) into 
aQ(k). This difference between QED on the one hand and YM theory and QCD on the 
other, precludes the use of unitary transformations to construct the set of physical states 
from the perturbative states {\n)} in the non-Abelian gauge theories; and that fact 
accounts for major differences between QED and the non-Abelian YM theory and QCD. 

In formulating a procedure for constructing the state {uo) in non-Abelian gauge theories, 
we note that whereas f2"(k)|z/) = and f2"*(k)|i/) = are two wholly independent conditions 
for QED — and the former alone suffices to define the Fock space for this Abelian theory 
— these two conditions are not independent for YM theory and QCD. In YM theory and 
QCD, Eq. (Q) requires that fl°-*{k)\h') = too The appropriate condition for imposing 

Gauss's law in YM theory therefore is not Eq. (H), but 

[n''{k) + n'*i-k)]\u) = o. (9) 

We will express this condition as 

[6^(k) + Jo«(k)]|z.)=0, (10) 

where we define 

6^(k) = fciK(k) + a-(-k)]. (11) 



We will transform the perturbative vacuum state |0) — the state that is annihilated by 
a°(k), ciQ(k) and ^^(k) — into a state |i/o) that implements Gauss's law. We represent |z/o) 
as a product of two operators acting on the perturbative vacuum |0), in the form 

Wo) = ^ H |0) , (12) 

where the operator product \E' H represents a non-unitary transformation operator. We also 
define a state |0o) = S |0), so that 

6^(k)|0o) = O. (13) 

Eq. ([13|) is satisfied by H = exp{— X^k; '2ij*('*)'2Q*(~'*)}5 as is confirmed by the observation 
that [aQ(k), S] = — (— k) S and [ag (k), S] = 0. The resulting state |0o) is not normaliz- 
able — it is essentially the "Fermi" vacuum state |]^, which is not commonly used in QED, 
but reappears here in the non-Abelian theory. 

The construction of \I' involves solving the equation {bqCk) + Jq (k)} \1/ \(f)o) = \l/6Q(k)|0o), 
or equivalently, 

rQik),^] = -JS{k)^ + B^Q, (14) 

where Bq represents any operator product with bq on its extreme right-hand side, so that 
^q\4'o) = 0- Eq. (|T^) is an operator differential equation, in which the commutator plays 
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the role of a generalized derivative. We introduce the following notation for the constituent 
parts from which ^ will be assembled: 

<(r) = A^v) , (15) 

^?(r) = AUr) , (16) 

and 

X'^ir) = [§At{r)] , (17) 

where, [af(r) +<(r)] = Af{r), and since yl£i(r) = ^^^^^"(r)] , xf{r) = diX'^{r). We also 
need to define the combination 

<).(r) = [af (r) + ^xf (r)] , (18) 

with T) integer- valued. 

We will also use the preceding operators to form the following composite operators: 

V'5).(r) = r^"' X'^ir) Qf,),(r) = A'-(r) [af (r) + i^f (r)] , (19) 

V'?2).(r) = -r'^' /'"^ A"^(r) Q(^).(r) X\r) = -/"/^^ /^^^ A'"(r) [af (r) + |xf (r)] A'^(r) , (20) 
and the general 77-th order term 

^J,).(r) = (-1)^-^ /Sf ^f.)(r) , (21) 

in which 

^f.)W= n A'-N(r), (22) 

m=l 

and 

jaf3y ^ yQ;[l]/36[l] j&[l]a[2]6[2] j6[2]a[3]6[3] . , , j&[??-2]a[7?-l]6[r)-l] jb[n-l]a[n]j ^23) 

where /('J^^ = /"^^. 

We also define the composite operator 

A^^igJdri;J,^,{r)U]{v), (24) 
which is useful because it has the important property that its commutator with ^^(k), 

[&^(k), A^] = -gr'^^J dr e-"^-- [af (r) + xf (r)] n7(r) 

-f /""^ J dr e-'^" A*" [ain7(r)] , (25) 
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generates — Jg (k) when it acts on the "Fermi" vacuum state |0o)- We observe that 

d,U]{r) = Y.bUk)e^^^-\ (26) 

k 

SO that diU] {r)\(pQ) = 0, and 

[6^(k), A]|0o) = -Jo(k)|0o). (27) 

We might expect that the simple choice \1/ = \l/o = exp(^i) would solve Eq. (|T^), but that 
expectation is not fulfilled. This is due to the fact that the commutator [^^(k), ^i] does 
not commute with Ai. The expression diU]{r), when it arises in the midst of an extended 
sequence of operator-valued factors, does not act on the Fermi vacuum state, and does not 
vanish. On its way to the extreme right of the expression, where it ultimately acts on 
|0o) and vanishes, diU]{r) produces extra terms as it commutes with ^I's; and that fact 
disqualifies \l/o ClS db solution of Eq. (|14D. 

To address the problems that arise because Ai and [^^(k), ^i] do not commute, we 
make the following modifications. First, we replace exp(^i) with || exp(^i) ||, where the 
\\0\\ designates a variety of 'normal order' in which all functionals of momenta, ^[Hj], 
appear to the right of all functionals of gauge fields, For example, in the n^^ term of 

II exp(^i) II the product || (A^)"^ \\ represents 

UAr\\ = i^9rl D(l,---,n)^-),^(l)^-)j2) ■■■^S,„(n) n-(l) n-(2) ■■■ n-(n) , 

(28) 

where D(l, ■ ■ ■ ,n) denotes dri ■ ■ ■ dtn, and the integer argument n in the ip^s and II's rep- 
resents r„. This 'normal order' has the effect that, in [^^(k), || exp(^i) ||], the diUj^r) 
produced by an integration by parts, appears among the II's and can annihilate the |0o) 
directly, moving only through other II's (or functionals of II's) with which it commutes. 
The {—g)f"''^'^J dr e~^^'^ [a'^{r) + x^(r)] needed as part of —J^ik) is created to the left of 
II exp(^i) II, as required; but the remaining n?(r) appears to the right of all the ip^s, with 
which it does not commute. As n?(r) moves to the left, to constitute the complete — Jg (k), it 
generates unwanted contributions that disqualify even || exp(^i) || as a solution of Eq. (|14|). 
To compensate for the failure of || exp(^i) || to satisfy Eq. (|T^, we extend Ai so that it is 
only the first term in the infinite operator- valued series A, given by ^ = J2^=i -^n, where 
the An with n > 1 will be given later in this section. \1/ is then given as 

^ = II exp(^) II . (29) 

All the An will consist of functionals of the gauge field, Af^r), multiplied by a single mo- 
mentum, n7(r), so that it becomes useful to express all the An as 

An = ^9''JdrAJ^),{v)U]{v). (30) 

Eq. (^) is useful as a definition of Aj^-^^^r), and as an identification of this quantity as a 
functional of gauge fields only — canonical momenta are never included in A7 ^^{r). 
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In Eq. ( [14D we can now replace \I/ with Eq. p9| ) to give 

[6^(k), II exp(^) II ] + Jo"(k)|| exp(^) || ^ , (31) 

where the symbol ~ is used to indicate that we have suppressed the state \4>o), that should 
appear on the right of all operator products. We can expand Eq. (|31|) into the form 

oo oo 

II [6^(k), exp(^X)] exp(^i) || + || [6^(k), exp(^i)] exp(5] A) || 

n=2 n=2 

+ Jo"(k)||exp(^)||^0, (32) 

where the || || -ordering eliminates further contributions from the Baker-Hausdorff-Campbell 
formula. Since [^^(k), || exp(^) || ] = || [&Q(k), A] exp(^) ||, we can use Eq. (|25D to eliminate 
jQ(k)|| exp(^) II and to rewrite Eq. (^) as 

oo „ 

II [fe^(k), ^ A]exp(^) II - \\gr^^ / rfre-^''-^[af(r) +xf(r)] [exp(^), n7(r)] || ^ . (33) 

n=2 

We can also use [|| exp(^) ||, n7(r)] = || [A, n7(r)] exp(^) ||, to giveQ 

oo „ 

II { [6^(k), E A] - 9r^' / rfre-'^-[af(r) +xf(r)] [A, UJiv,)] } exp(^) || ^ . (34) 

n=2 

After simplifying Eq. (|3^) and expanding, we establish 

oo „ oo 

[&Q(k), E A] - gr'-' / dre-^'^-^[aUr)+xUr)] E A, n7(r)] ^ , (35) 

^ O 1 



n=2 n=l 



as a sufficient condition for the validity of Eq. (0). We now rewrite the limit of the sum in 
the second term, J2'^^=iAn =^ J2^=2-^n-i, to give 

oo „ oo 

[&Q(k), E^n] - 9r^' / rfre-'^-[af(rO +xf(r)] [E A-i, n7(r)] ^ . (36) 

n=2 n=2 

Requiring Eq. (^) to hold for all values of g, we obtain the the recursion relation 

[6^(k), An] ^ g r^'J rfr e-^*^-^ [af (r) + xf (r)] [X-i, n7(r)] , (37) 
which holds for n > 1, because = 0. 

We have been able to construct the first six terms of the A series and to confirm their 
consistency with Eq. (P7|). Because of the structural regularity of Ai - Aq, we can also infer 
the form of the general An- The expressions for these An are most easily given in a partially 



^in an expression || [i^,C]CII) the commutator is always to be evaluated before the double-bar 
ordering is imposed. 
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recursive way, in terms of the ^^^^^(r) previously defined in Eq. (^). The definition of each 
An (with n > 1) contains references to Aj^,^j^{r) with n' < n, and, in turn, together with 
Eq. (pOD, defines the new ^^^^.(r). The terms in the A series are given by Eq. (0) and by 

A, = '4 j dv i;l^^{r) n/(r) 
+ ig' r^^Jdr |[-4^,),(r)] af (r) n7(r) , (38) 

^3 = ¥/rfr^J),(r) n7(r) 
+ ^9'r'''Jdr^[Af,^^{r)] af(r) n7(r) 

+ zg' r^'<j dr |[^fi),(r)] [6,, - if^]^f,),(r) n7(r) , (39) 

A, = 4jdv i^l^^ir) n7(r) 

+ ^^?V"^Y^r|[^f3)^.(r)]af(r)n7(r) 

+ zg' r^^J dr |[^f,),(r)] [5,, - i^] ^f,),(r) n7(r) 

+ zg' r^^J dr ||[^f,),(r)] [5,, - i^] ^f,),(r) n7(r) 

+ 4 /''V «f (r) n7(r) , (40) 

A, = '4 j dr i;l)Ar) n7(r) 

+ ^/r^Y^r|[^f4)^.(r)]af(r)n7(r) 
+ r^-'j dr |[^f3),(r)] [5.. - 1^] ^fi),(r) n7(r) 
+ r^-'j dr |[^f2),(r)] [5.. - \^] A%^,{r) n7(r) 
+ r^-'j dr |[^fy,(r)] [5,, - i^] ^f3),(r) n7(r) 
+ f |[^f2),(r)]||[^Ji),(r)] af (r) n7(r) 

+ 4 |[^(i),(r)]|[^f2).(r)] af (r) n7(r) 

+ 4r''f''Jd^ [^^^ - n7(r) , (41) 

and 

A = f/rfrV^J),(r) n7(r) 

+ r^Y^r |[^f5)^.(r)] af (r)n7(r) 
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+ tg' r^^j dr J[^^3),(r)] {5^k - ^f2),.(r) n7(r) 

+ ^g' r^^j dr |[^^,),(r)] [5,, - if^] A^^.^.ir) n7(r) 
+ zg' r^^J dr |[^^i),(r)] [6,, - if-] ^f,),(r) n7(r) 

+ f /''Y &[-^(3b (r)]i['^(i)'^(r)] «f (r) nZ(r) 

+ ¥ r''^ f'^J |[^f2),(r)]i[^f2).(r)] af (r) n7(r) 

+ ¥ r'"' /''Y |[^(i),(r)]|l[^?3).(r)] af (r) n7(r) 

+ '4r''f''jd^ |[^{2),(r)] il[^fi),(r)] [5. - if^]^fi),(r) nZ(r) 

+ '^r''f''Jdr |[^fi),.(r)] [5. - if^]^fi)Kr) n7(r) 

+ '^r''f''Jdr |[^fi),(r)] il[^Ji),(r)] [5, - l^]A'^,),{r) U]{r) 

+ f r'^' f"" r'^J dr |[^^,),(r)] ||[^fi),(r)] §[Aly{r)] af (r) n7(r) . (42) 

To arrive at a form for An for arbitrary n, it is convenient to define A]{r) by Eq. (|30D 
and by 

oo 

Mr) = T.9''AJ^)^i^), (43) 

n=l 

SO that 

A = tJdr:^{r)U]{r) . (44) 

We also define 

B^,{r) = {af (r) + [6,, - ^¥^]:4f(r) } , (45) 

and 

^U^) - n = |[^(r)]|[:4?^(r)]---§[Xf^(r)]||[:4^^ (46) 

m=l 

where t] represents the same integer- valued parameter originally introduced in Eq. (|18|), 
that we will now observe to be necessary for matching each B^^y {r) with its corresponding 
A4'^^^{r) in the integral equation for A]{r) given below. We now formulate this non-linear 
integral equation for Af{r) as 

oo „ 

^=T.f drrM^i^) U]{r) 

ri=l 



9 



+ig fft; j dvMUv) W^^ir) n7(r) 
+ ^4 ff^^ j dv M%){v) B^^ir) n7(r) 



/^rA<f,)(r)i3f^),(r)n7(r) 
+ ••• , (47) 



or, more succinctly, 



=1 



-4 = E^/^r{V^;;jr) + drMUv)BUr)}U]ir). (48) 



We observe that the leading terms of the perturbative solution of Eq. ( ^71 ) — or equivalently 
Eq. p8|) — agree with Ai - Aq, the explicit forms given in Eqs. (0) and (|38|)-(p2D. We have 
confirmed that the entire perturbative series — An for arbitrary n — correctly satisfy the 
recursion relation given in Eq. (p7D; but we have not yet established that fact with complete 
rigor in so far as concerns its extension beyond n = 6. Previously in this paper, we have 
shown Eq. (|37|) to be a sufficient condition for the implementation of Gauss's law. 



|z/o) is not the only state that implements Gauss's law. Any state |z/k) = ^ '=''^sUk)|0) or 
|^'ki --ki) = ^1' H a°^^'(ki) • • • a"^''''(kj)|0), where a"^(k) creates a transversely polarized gluon, 
is annihilated by the Gauss's law operator Q. The important question then arises: What 
changes occur in the S-matrix when the states |z/ki- -ki) are substituted for the perturbative 
states |nki...ki) as incident and scattered states? We will not discuss this question in detail 
in this paper. But we observe that, while it has been shown that there is no change in 
S-matrix elements when the |z/) and the corresponding |?t,) states that they replace are 
unitarily equivalent [^0, in this case, in which the transformation is not unitary, there 
willhe changes in the S-matrix elements when the \i>) states are substituted for the \n) states. 
These changes, however, do not appear in the lowest order of perturbation theory, since the 
non-unitarity does not arise in the leading term of \E'. The most interesting possibility, of 
course is that the contribution of as a whole can be evaluated, and the non-perturbative 
effect of Gauss's law on the S-matrix assessed. 



Another question to be addressed deals with the non-perturbative solutions of the non- 
linear Eq. (^Tf ) or equivalently Eq. (|48|). Further work is required to clarify how these 
non-perturbative solutions are related to the gauge sectors connected by the large gauge 
transformations [§]. 

This research was supported by the Department of Energy under Grant No. DE-FG02- 
92ER40716.00. 
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